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Analytical solutions are developed for tracing radio rays through a spherically stratified iono- 
spheric layer. The electron density profile is represented by a series of parabolic sections, each sec- 
tion being bounded by levels of maxima or minima of electron density and by the middle values between 
these extremes. 

With this model, relatively simple functions are obtained which give the length of the radius vector 
in terms of these levels. For example, it is shown that, under certain conditions, the ray path is a 
simple ellipse or a quasi-ellipse. The cases of reflection and penetration are considered for waves 
originating on the earth's surface and for waves incident from outside the ionosphere. 



1. Introduction 

The need for an analytical formula capable of ex- 
pressing the path of a radio wave in the ionosphere was 
raised at the meeting on Oblique Sounding of the 
Ionosphere [IQSY Notes, 1963] held at Lindau, Ger- 
many, in May 1963. This led to the derivation of 
equations which express the following kinds of radio 
wave propagation: (a) penetration of the layer, (b) 
oscillation along a circular path, and (c) reflection. 
Although the formulas refer to propagation in the ab- 
sence of the earth's magnetic field, they may be ap- 
plied also if a slowly varying magnetic field is present. 



2. Propagation in Perfectly Spherical Ionized 
Layers 

A "perfectly spherical" layer is assumed to be one 
in which the levels of constant electron density (N) 
are concentric spheres. Each sphere is defined by 
its radius of curvature (r) (fig. la). The variation of 
electron density from level to level is quite arbitrary; 
it must, however, remain the same in the whole region 
investigated [Woyk, 1964]. 

For waves propagated in spherically ionized layers, 
there are two fundamental formulas defining the ray 
path: Snell's Law, 



r n sin i = constant = r m sin i c — r t — n r 00 



r, n r = rK sin ie, 



(i) 




FIGURE 1. Geometric representation of symbols used in the text. 
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which is valid for any point on the ray path and the 
refraction formula 



f „ dn tan i , 

dK = ; dr, 

dr n 



(2) 



which expresses the deviation of a ray from its pre- 
vious direction. The meaning of the symbols is best 
seen from figure la. (r t ) is the sphere to which the 
ray path is tangent whenever it is outside the ionized 
layer (iV=0, n=l), while (r/ n) ) is the radius of the 
sphere to which the ray is tangent whenever it is within 
the ionized layer. 



n 



(3) 



where (n) is the phase refractive index 



n = Vl- 






r-$. n. 



:N n 



(4) 



where (/) is the wave frequency and (f c ) is the critical 
or penetration frequency of the layer. 

Three categories of radio propagation must be taken 
into account [Woyk, 1963]: 

(a) waves penetrating the layer, [Woyk, 1958 and 
1963] 

(b) waves propagated along circular paths, 
[Chvojkova, 1954] 

(c) reflected waves [Woyk, 1959] . 

(a) For waves penetrating the ionized layer, the fol- 
lowing useful theorem has been derived [Chvojkova, 
1958 and Woyk, 1963]. Both linear parts of a ray 
path — before entering the layer and after having left 
it — are merely turned about the center of curvature of 
the layer by the angle (/?), which can be expressed in 
the form 



R = 



sin i c 



fNdr 



2 cos 3 i c N m r m 



(5) 



where (/?) is in radians. From this equation we see 
that angle of refraction of the ray penetrating a layer 
does not depend upon the shape of electron density 
profile, but only upon the total electron content in a 
vertical column of unit area. 

This formula is accurate to (q) decimal places only if 



cos* i c 



10-« and A ^ r m 10"^, 



(6) 



where A is the thickness of the layer. 

(b) Radio waves can be propagated along two cir- 
cular paths with radii r x and r' x (fig. Ic), defined by 



dn 
dr' 



(7) 



The upper one (r^) lies just below the maximum of the 
layer (r m ). Because of a strong dissipation of the beam 
propagated along level (r^), this type of (high ray) 
propagation is not of much practical use. The propa- 
gation along the lower one of radius (r'J lies just above 
lower boundary of the layer and is very stable. Rays 
which deviate slightly from this circular path are 
forced to return back to it perpetually either by re- 
flection from the denser upper regions or because of 
the curvature of the layer. Thus the wave propagates 
alone r' x in small oscillations. 

(c) Formulas derived for reflected waves [Woyk, 
1959] are either so complicated or so simplified that 
they are of no practical use. Therefore it has been 
necessary to search for a suitable electron-density 
distribution which would approximate the real one and 
which when introduced into the refraction formula 
would lead to an analytical solution. 

3. Analytical Ray Path Formulas 

Let us divide the ionosphere into zones as shown 
in figure 2. Each zone is bounded by spheres of radius 
r x or r s . The r x 's denote levels of local maxima or 
minima of electron density, more specifically denoted 
as tm and r^, respectively, and r s lies in the middle 
between the two extremes, such that 



rs=h(r M + riL). 



(8) 



The actual electron density at a local maxima is given 
by M and at a local minima by /x, with X representing 
the electron density at either extreme and that in the 
middle given by 



S=i(M + /i). 



(9) 



The maximum electron density of the entire ionosphere 
is further specified as N m and r m refers only to this 
level. 
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Figure 2. Division of an ionized layer into appropriate zones. 
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The shape of the electron density profile of each 
zone can be approximated by a quasi-parabola having 
this form: 



To obtain an expression in terms of r and <p, dR and 
di must be introduced in (19) as functions of r only. 
From SnelFs Law we have 



N=X+(S-X)\ 



r* r-r x \ 2 
- — r.v r I 



(10) 



The whole shape of the ionosphere can be repre- 
sented by a series of such parabolas which all together 
give a smooth curve with the exception of the junctions 
at the levels r s . 

If we introduce (10) into (4) and make the following 
substitutions, 



W- 



F(XS)/ r s 



N m -FX\r s -r Xi 
n 



V 



1-™ 

N„, 



(11) 

(12) 



we have 



n 2 = \ + W 



r., 



{nrf = r\W+ 1) ~ Wr x r + Wr 2 
Consequently, it is also true that 

n = 



rt 



v 



(13) 



(14) 



(15) 



N m 



With a simple nocturnal ionosphere (fig. 2b), it is 
appropriate to use 



N = 



-N n 



U + 1)- 



r—r x 



r x - 



and therefore 



w=+— 



2 + {\ + l)F\r s -r x 



r s 



(16) 



(17) 



where the upper (minus) signs refer to the outer (con- 
cave) zones I and V and the lower (plus) signs refer to 
the inner (convex) zone II and IF. 

Since we wish to derive an analytical formula, it has 
been found appropriate to introduce polar coordinates 
r and (p. It follows from figure lb that 



and 



<p = R — (i2 — i\) 
d(p = dR — di. 



(18) 



and 

Ji 2 r 2 cos 2 i = (nr) 2 —Tf. 

By differentiating (13), we get 

dn = Wr x {r — r x ) 
dr nr 3 

and substitution in (2), 

,r> Wr x (r-r x ) 

dK = — r t - dr. 

r(nr) 2 \{nr) 2 - 7f 

Rearranging (20) and differentiating, we have 

,. di. . _ (W.+ \)r-Wr x y 

di= — di=—r t ; : dr. 

dr (nr) 2 V(nr) 2 -r 2 

Introducing (23) and (24) into (19) we get 

dr 



d(p = n 



rVrW+l)-2Wr x r+Wrl-rf 



(20) 



(21) 



(22) 



(23) 



(24) 



(25) 



which has an analytical solution. Depending upon 
initial conditions, there are three solutions that must 
be distinguished. 

Case I. For W <\ — I , which is fulfilled in concave 

w 

zones I and I', the solution of (25) is 



1 -he sin k((f — (po) 



(26) 



where 



Wr 2 -r 2 V(JF+ l)r 2 -rr 2 , / 



Wri 



(27) 



Since the reference for the angle <p has not as yet been 
chosen, let us set <po = — 90° (at reflection <p will be 
zero). Equation (26) can now be written as 



r= 



1 — e cos k(f 



(28) 



If 4=1, the equation (28) is that of a simple ellipse. 
Usually A: is a large value and thus (28) is a quasi- 
ellipse whose (p has been A:-times contracted as indi- 
cated in figure 3a. For € = the path is circular and 
the wave propagated along r x . 



(19) Case II. For r> 



£f. 



which is fulfilled in the con- 
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vex zones II and II' two other circumstances must be 
distinguished. 

If e 2 < 0, the solution of (25) has the form 



1 — e' sinh k'<p 



where 



e' = \T^ and k' = V~k 2 . 



(29) 



(30) 



In this instance the ray passes the local layer without 
reflection as shown in figure 3b. 
If e 2 > 0, the solution becomes 



1 — (±e) cosh k'ip 



(31) 



as represented in figure 3c, and the wave is reflected. 
If we take e positive, we are dealing with a wave com- 
ing from above the ionosphere and reflected outwards. 
A negative e corresponds to waves from below reflected 
downwards. For e = 0, the wave propagates along 
the circular path r x . 

Any ray path is sufficiently determined by F, r x , 
r s , and r t . From these parameters one can find the 
values of p, e, and A:, and hence one can determine 
all the partial curves of which the whole path is com- 
posed. Although composed of several fragments, 
the trajectory should be a smooth curve with a pos- 
sible exception at r s . 




The value of r t must often be deduced from any one 
of several considerations: the angle of incidence at 
the earth's surface, the level of reflection, or the 
angle at a level r, or by the sphere r* n) to which the 
ray is tangential whenever it passes a level with the 
index of refraction n. In any of these cases, r t is 
easily determined from SnelFs Law. 

4. Some Additional Remarks 

The above formulas can also be used when a mag- 
netic field is present. If the space variation of the 
field is not great, it is sufficient to introduce the 
Appleton-Hartree formula into W . 

The level of reflection r r is found from (28) and (31) 
by introducing cp = 0. Thus 



r r 



= P 



(32) 



No ray reflected below r^ is reflected down. All are 
deviated into the denser upper regions, since they are 
less curved than the earth. 

The circular paths are defined by e = 0. The r x or 
r^ are defined by r — p or more specifically, 



roc or r 



"r+i 



(33) 



It follows that one circular path (r^) lies in zones I 
only, just above the minimum electron density or lower 
boundary of the layer. The stability of this type of 
propagation is great as is evident from figures lc and 
3a. A deviation from the circular path brings the ray 
into oscillations along r^, as previously mentioned. 
In zones I', no circular paths can occur and all rays 
reflected there are reflected outward. 

The other circular path lies only in zone II, closely 
below the maximum of electron density. Figures 
lc and 3c show that these paths are unstable and the 
ray is easily dissipated on both sides. 

With increasing /, \W\ decreases. r x and r^ ap- 
proach each other slowly, drifting slightly toward 

\r s — r x \ 
r s . For thin layers, — = 10~ 3 , they come together 



al 



30. 

/ 

fr 



For thicker layers, 
= 5 (see table 1). 



= 2.10" 2 , they 



FIGURE 3. Representation of ray paths for the analytical solutions. 



tt-=30. For thicker layers, — ~ 

Jc r x 

join at 

When a sharp gradient of electron density is to be 
investigated, formulas (28) and (31) or (29) can be used. 

ror expressing -7— ~ °°, however, we must set r x — r s 
at 

~ 0, or M and /x much above N m and much below A^ = 0, 

respectively. This is consistent since the parts above 

N m and below A^ = are not investigated. When 

-7— increases very much, \W\ and \k 2 \ also increase 

very much. The oscillations of figure 3a become more 
squeezed and the ray path is much sharper at the 
level of reflection. 
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Table 1 



Mr 


3:1 


5:1 


7:1 


10:1 


100:1 


r s fr m 


W in Zone I 


0.999 
.990 
.980 


-55,400 
-545 
-133 


-19,960 
-196 
-48 


10,200 
-100 

-24.5 


-4,990 
-49 
-12 


-49.9 

-0.5 
-.12 


r s /r m 


W in Zone II 


0.999 
.990 
.980 


62,400 
612 

150 


20,800 
204 
50 


10,400 
102 

25 


5,050 
49.5 
12.1 


50 

0.49 

.12 


rjr,„ 


W in Zone II' 


1.001 
1.010 
1.020 


62,500 
638 
162 


20,800 
212 
54 


10,100 
106 

27 


5,050 
51.5 
13.1 


50 

0.51 

.13 


r s /r m 


W in Zone I' 


1.001 
1.010 
1.020 


-55,600 
-567 
-144 


-20,000 

-204 

-52 


-10,200 
-104 
-26.5 


-5,010 
-51 
-13 


-50 
-0.51 
-.13 


rjr m 


r'*lr,„ in Zone I 


0.999 


0.998 
.982 
.967 


0.998 
.985 
.980 


0.9981 
.990 


0.9982 




.990 




.980 














r s lr tn 


r^lr,„ in Zone II 


0.999 


1.0000 

0.998 

.993 


1.0000 

0.995 

.980 


0.9999 
.990 


0.9998 




.990 




.980 














r r 


W in Zones II or II' 





0.0625 


0.0208 


0.0104 


0.00505 


0.00005 


r T 


W in Zones I or I' 





-0.0555 


-0.0200 


-0.0102 


-0.00500 


-0.00005 



Some apparent discrepancies arising during an 
analysis are usually avoided if we consider that we 
wish to investigate only the regions bounded by r x 
and r>. Outside this region, the theoretical curves 
increase rapidly above N m or below N = Q, and the 
region I goes over into I' at the level r x . Thus out- 
side the boundaries r x and r s , the result no longer 
represents the investigated case. 

At W~0 attention must be paid to the fact that in 
the zones II and II' (28) is valid even for positive 

W< (— V while (29) and (31) are true only forJF > (y- 



This condition exists only when f> f c or when one of 
the boundary levels approaches r x = 0. Since this 
eventuality is of no practical importance for terres- 
trial radio communication it will be analyzed else- 
where [Woyk, 1965]. 

It is surprising that in a zone bounded by r x = 0, the 
ray path is rectilinear. This is because (10) gives a 

constant electron density N = —z~ in the entire zone, 

except where r = 0. 

For the sake of a rapid survey, some values of W 
and r x have been given in table 1. They correspond 
to a simple layer similar to figure 2b but symmetrical 
about r m . For r# = 0, the values of W refer to (k); 
they are independent of r s . When r x ~ 0, (7) corre- 
sponds to a layer much different from that in figure 2b. 
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this paper. Further I thank Dr. K. Davies, Dr. R. G. 
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